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Abstract 



For one-dimensional phase-turbulent solutions of the Kuramoto- 
Sivashinsky equation with rigid boundary conditions, we show that there is a 
substantial variation of the correlation time t c (x) with spatial position x in 
moderately large systems of size L. These results suggest that some time- 
averaged properties of spatiotemporal chaos do not become homogeneous 
away from boundaries for large systems and for long times. 
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For spatiotemporal chaotic systems (whose spatial correlation length £ is smaller than 
the system size L an interesting question arises as to whether time-averaged properties 
may vary spatially. In a thermodynamic limit of infinite system size and for averages over 
an infinite time interval, one would expect time-averages to be spatially-independent as a 
consequence of ergodicity and of rotational and translational invariances of the underlying 
equations. For experiments or simulations of finite-size, however, spatial variation may arise 
from boundary conditions, from broken symmetries, or from a combination of both. 

Time-averaged inhomogeneous patterns have indeed been found for finite-sized systems 
both numerically and experimentally. Zaleski [[J and Pumir studied long-time statistical 
properties of spatiotemporal chaotic solutions of the one-dimensional Kuramoto-Sivashinsky 
(KS) equation 

v* + u xx + u xxxx + uu x = 0, for < x < L, (1) 
with rigid boundary conditions 

u = u x = 0. (2) 

The time-averaged pattern p(x) = (u(x,t)) [F2] and simple statistical quantities like the 
kurtosis |§ (derived from time series measured at a single point in space) were found to vary 
spatially, but only for spatial intervals small compared to the system size and only near the 
boundaries. An example is given in Fig. [I], which represents a more thorough time-averaging 
than that reported in Ref . . We found that the maximum magnitude of the spatial average 
was independent of the system size L over the range 50 < L < 2000, suggesting that the 
mean spatial pattern is a consequence of the boundary conditions. 

More recently, Gluckman et al 0] found quite striking time-averaged two-dimensional 
spatial patterns in a Faraday crispation experiment for system sizes that were about five 
times larger than the correlation length £. Unlike the results found for the KS equation, 
the time-averaged experimental patterns showed a highly-ordered periodic lattice with the 
symmetry of the lateral boundaries (square or circular). Similar results have been reported 
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by Ning et al for rotating convecting fluid in cells of different symmetry ||. Although a 
systematic study of how the magnitude of these patterns decreases with aspect ratio was 
not made, the results for both numerics and experiments are consistent with the expected 
behavior that the magnitude of the time-averaged patterns should vanish in the infinite- 
aspect-ratio limit. 

Somewhat related to these investigations of time-averaged patterns are two studies |||7j 
concerning whether the fractal dimension D 2 [§ varies in space. This is a possibility since 
the fractal dimension is related to the strength and number of active modes, and these can 
vary spatially for finite-aspect-ratio systems 0. Sato et al || and Ciliberto obtained their 
results for chaotic convection in rectangular cells of modest-sized aspect ratios, respectively 
15 x 1 and 4x1. They measured long time series at many spatial points and analyzed 
these time series with the Grassberger-Procaccia algorithm [10] to extract D 2 (x). Sato et al 



found a strong spatial dependence of Fourier modes but no clear spatial dependence of D 2 . 
Ciliberto found a stronger spatial dependence of D 2 in the smaller cell. Both experimental 
results were not entirely convincing since dimension algorithms based on time series are 
difficult to apply in a reliable way [II]. Although an unambiguous dependence of fractal 



dimension on position was not found, Ciliberto did observe a strong spatial variation of the 
root-mean-square of an observable, with the variation occurring on lengthscales smaller than 
a roll size. 

We note that the above studies concern issues that are somewhat different than those 



raised by the study of coherent structures in more fully developed turbulent flows ||i2| , |T3 
Coherent structures typically persist only over finite time intervals and long-time averages 
have not yet been carried out to see how the time-averaged pattern may be related to existing 
coherent structures. 

In this Letter, we look at a different aspect of spatial-variations of time-averages. We 
investigate the question of whether temporal correlations, as quantified by a correlation 
time r c , can vary spatially. This is a rather different question than looking at time-averages 
of fluctuations since time correlation functions involve comparing information separated at 



time intervals that can be quite large. In the following, we show empirically that there 
is an unexpected spatial variation of the correlation time for a one-dimensional model of 
spatiotemporal chaos but only for rigid boundary conditions and only for phase turbulent 
states |b|]. Unlike the numerical |2|,[|] and experimental |||| results discussed above, the 
spatial dependence of r(x) is strongest in the interior, away from the boundaries. This 
spatial dependence exists only in the phase-turbulent KS equation, is weakly sensitive to 
the choice of initial conditions, and disappears with sufficiently long averaging in the more 
disordered defect-turbulent regime of the CGL equation [|l5| , p^1 . 

One implication of these results is that not all features of large-aspect-ratio spatiotem- 
poral chaos become homogeneous away from boundaries, even for systems that are quite 
large compared to the correlation length (for our calculations, this ratio is about 20). Ex- 
perimentalists and computational scientists will have to be wary of which point they pick in 
space to sample information, at least in systems with long-lived temporal correlations as is 
the case for phase turbulence. 

We obtained our results by numerical integrations over times T as large as 10 6 time units 
for Eqs. (|I|) and (Q), and for the closely related but more general complex-Ginzburg-Landau 
(CGL) equation M 



d t A{x,t) =A + (l + i Cl )d 2 x A-(l-ic 3 )\A\ 2 A, (3) 

with Dirichlet boundary conditions 

A(0, t) = A(L, t) = 0. (4) 

The CGL field A(x, t) is complex- valued and the parameters c\ and C3 are real- valued. The 
KS equation can be derived from the CGL equation in the Benjamin- Feir unstable 
limit C1C3 — > 1. For fixed C\ > 1.8, there is an apparent transition from phase turbulence, for 
which the field magnitude \A\ is bounded away from zero for all time, to defect turbulence, 
in which isolated space-time defects of zero field A value occur [TjJ. It is not yet known 



whether the two phases are distinct in the thermodynamic limit of an infinitely large system, 



e.g., phase turbulence may simply be a defect turbulence with an unobservably low but finite 
density of space-time defects. Unpublished results of ours, extending Shraiman et al to much 
larger system sizes, suggest that the phases are in fact distinct. 

We integrated Eqs. (JTJ) and (0) and Eqs. (|3|) and (f|) with second-order-accurate finite- 
difference codes with time-splitting of the linear and nonlinear operators |T7|] . At each of 
many coordinates x on a regular spatial mesh spanning the interval [0, L], we calculated a 
correlation time t c (x) from the x-dependent time correlation function 

C(t; x) = ^ (u(t + r, x) — u(x)j (u(t, x) — u(xfj J , (5) 

where u(x) is the function representing the time averaged mean at each point in space and 
where the brackets indicate time averaging over information at position x. This function was 
found to decay approximately exponentially for small times. We used the root-mean-square 
width of C (r; x) up to the first zero crossing to determine the correlation time r. Our results 
(especially Fig. Rf) are not sensitive to our method of extracting r, e.g., we found similar 



results when estimating r by the first-zero-crossing of C(r) or by the integral time scale [|T8 

We note that the long-time long-wavelength dynamics of the KS-equation and of the 
phase-turbulent regime of the CGL equation is believed to be described by the Kardar-Parisi- 
Zhang (KPZ) equation [|19H 22|1 which describes the roughening transition of an interface. If 



this asymptotic description is correct, the time correlation function in one spatial dimension 
does not decay exponentially with a single well-defined time scale, but instead decays as 
a stretched exponential of the form exp(— Bt 2/3 ) with [3 = 1/3. We ignore this subtlety in 
what follows since, as shown in Fig. 0, there is a substantial decay of correlations over a 
finite time interval and the system sizes are small compared to the lengthscale over which 



the KPZ description becomes valid pOfl 



We also ignore a possibly serious issue concerning whether our numerical results are sta- 
tistically stationary so that all transients have decayed sufficiently Naively, one would hope 
that an integration time T of order 10 6 that is many orders of magnitude greater than the 
correlation time r c < 40 would be sufficient for transients to decay. However, some calcula- 



tions by Shraiman |23| suggest that KS phase turbulence with periodic boundary conditions 



is transient, with a decay time that grows exponentially with system size L. (It is not known 
if this conclusion generalizes to the boundary conditions Eq. (|2|) although numerical calcula- 
tions suggest that this is unlikely). The possibility of extremely long transient times has also 
been argued from numerical studies of coupled map lattices pi| , |2"5| . We do not know of any 
quantitative method for distinguishing such long-lived transients from nontransient chaos, 
nor do we know whether there are observable dynamical implications of averaging over such 
a transient rather than over an attractor. We present the results below with the hope that, 
if they are transient, the slow decay times lead to an ergodic average that is approximately 
correct. 

Our results are summarized in Figs. |2| through [5|. Fig. ^| shows parts of the temporal 
and spatial correlation functions averaged over all space and time, from which correlation 
times and lengths were calculated for the KS equation with rigid boundary conditions. We 
estimate a value of r = 2 units and of £ = 7 units respectively and choose integration times 
and system sizes much larger than these values in what follows. 

Fig. |3|(a) shows the central result of this Letter. The correlation time of the KS equa- 
tion, as measured by the root-mean-square width of the temporal correlation function for 
rigid boundary conditions, varies substantially with spatial position. The magnitude of 
this spatial variation is quite large, with some peaks being as large as five times the average 
background value of about 8. The peaks show structure that varies substantially over length- 
scales comparable to the correlation length. Fig. |3](b) shows how the temporal correlation 
function itself varies over a small spatial interval 25 < x < 30 of the system. The variation 
in correlation times between x = 25 and x = 30 is caused by the presence or absence of 
the small dip in the correlation function around r = 12. The variations in the correlation 
function at large r remain about the same magnitude out to the largest values of r that we 
tested. When averaged over space, these variations disappear as shown in Fig. 0(a). 

From similar numerical calculations, we have found empirically that the positions of the 
peaks in Fig. ||](a) do not change with the length of the integration time T, with different 
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initial conditions or for changes in the numerical time resolution. However, the relative 
heights and widths of the peaks do change under these conditions. 

To obtain more insight about the origin of these peaks, we have made similar calculations 
for the one-dimensional complex Ginzburg-Landau equation, Eq. (|3]), which has a phase- 
turbulent regime similar to the KS equation but also a more strongly disordered defect- 
turbulent regime for which temporal correlations decay much more rapidly [|l^,^]. The 
spatial variation of Fig. |||(a) is not seen for defect-turbulent states with rigid boundary 
conditions, Eq. (|J), as shown in Fig. f|. Fig. f|(a) shows the time averaged mean of the 
amplitude of the field A as a function of position while Fig. |](b) shows how the correlation 
time varies with position. In both cases, peaks are observed near the boundaries while the 
bulk variations are small and smooth. The asymmetry in the heights of the outermost large 
peaks may be a consequence of some small initial asymmetry in the initial conditions. (We 
have not studied rigid boundary condition solutions in the phase-turbulent regime of the 
CGL equation, which presumably would give results similar to Fig. |3|.) We deduce that the 
spatial variation of the correlation time is a consequence of the long-lived time correlations 
present in phase turbulence. 

We have also tested the importance of rigid boundary conditions by numerically integrat- 
ing the Id CGL equation with periodic boundary conditions in the phase-turbulent regime. 
In this case, we find that there is no spatial variation of the correlation time as shown in 
Fig. |5|. Essentially identical results are found for the KS equation with periodic boundary 
conditions with the difference that one finds substantial spatial variations for short integra- 
tion times that slowly decay away in magnitude for longer integration times. These results 
confirm our earlier summary of previous numerical and laboratory experiments that 
average spatial patterns are driven by the presence of boundary conditions that break the 
translational symmetry. 

In conclusion, we have looked at a variation of the idea that time-average spatiotemporal 
chaotic patterns can have interesting mean structure. By numerical simulations on the KS 
and CGL equations in one-space dimension, we have found that the correlation time r c can 
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vary substantially in space when there are long-lived temporal correlations (phase turbu- 
lence) and for boundary conditions that break translational invariance. For dynamical states 
with more rapidly decaying correlation functions (defect turbulence) or for periodic bound- 
ary conditions, there is no interesting spatial variation of r c provided that one integrates 
sufficiently long in time. It would be interesting to explore these issue further in laboratory 
experiments on fluids and flames. 
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FIGURES 

FIG. 1. Spatial dependence of (a) the time-averaged mean-value (u) and of (b) 
the time-averaged square-deviation ((u — {u)) 2 ) for solutions of the one-dimensional Ku- 
ramoto-Sivashinsky equation Eq. (||) with rigid boundary conditions Eq. (||) on an interval of 
size L = 100. The temporal resolution was At = 0.1, the spatial resolution was Ax = 0.5, and 
the total integration time was T = 10 6 time units. For initial data, we used uniformly distributed 
random numbers in the interval [—0.03,0.03]. 

FIG. 2. (a) Plot of the temporal correlation function Eq. (||) averaged over all space 
and time for a numerical solution of the KS equation, Eqs. (||) and ((H). There is a sub- 
stantial decay of correlation over a time scale r c = 2. The numerical parameters are 
the same as those given in Fig. |]. (b) Space- and time-averaged spatial correlation func- 
tion S(y) = (u(t, x + y) — u(t)J (u(t, x) — u(tfj ^ ^ plotted over the interval < x < 50 for 
a system of size L = 100. There is an oscillatory exponential decay of correlations, defining a 
spatial correlation length £ = 7 that is small compared to the system size L. 

FIG. 3. (a) Plot of the root-mean-square (rms) width of the correlation function C(r) as 
a function of position x, for the same parameters used in Figs |l] and ||. A substantial spatial 
dependence of r c is found, (b) Details of the variations of C(r) over the spatial interval 25 < x < 30. 
The initial decay of the correlation function has not been plotted to allow smaller features of the 
function to be seen on the same scale. 

FIG. 4. (a) The time-averaged mean value and (b) the rms width of the time correlation 
function C(r) of the amplitude of the field in the CGL equation with rigid boundary conditions, 
Eqs. (0) and (Q). Results are presented for the parameter values c\ =3.5 and C3 = 1.5, which is 
in the defect-turbulent regime | |16| |. Spatial variations are found only near the boundaries. The 
temporal resolution was At = 0.05, the spatial resolution was Ax = 0.25, and the total integration 
time was T = 4 x 10 6 time units. 
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FIG. 5. Plot of the rms width of the time correlation function C(r) for periodic boundary- 
conditions of the Id CGL equation for parameter values (cl,c3) = (3.5,0.7). The numerical 
parameters were L = 100, 256 Fourier modes, At = 0.2, and T = 10 6 time units. After a 
sufficiently long time-averaging, the correlation time is no longer varying spatially. 
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